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Abstract 

We present the coefficient of the chromomagnetic interaction opera¬ 
tor, the only unknown coefficient in the Heavy Quark Effective Theory 
(HQET) lagrangian up to the 1/m level, with the two-loop accuracy by 
matching scattering amplitudes of an on-shell heavy quark in an exter¬ 
nal field in full QCD and HQET, and obtain the two-loop anomalous 
dimension of this operator in HQET. 


1 Introduction 

The leading order HQET lagrangian [1] 

Lq = Q v ivDQ v (1) 

has a unit coefficient by construction (here Q v = jQ v is a static quark field with 
velocity v). At the 1/m level two new terms appear [^]. || 

L = L 0 + Q v ((vDf - D 2 ) Q v + 9 ?Mq G a ^ Qv ( 2 ) 
2 m v 4?n 

(where the composite operators are normalized at the scale /x). For the coeffi¬ 
cient of the kinetic-energy operator 

Ckip) = 1 (3) 

holds in all orders of perturbation theory, due to reparametrization invariance of 
HQET B10- Only the coefficient of the chromomagnetic interaction 

operator is not known exactly. It can be found by matching scattering ampli¬ 
tudes of an on-shell heavy quark in an external chromomagnetic field in QCD 
and HQET up to 1/m terms. This was done in Q at the one-loop level; the 
one-loop anomalous dimension of the chromomagnetic interaction operator is 


1 





therefore known [^]. It is natural to perform matching at g « to, where C m (/i) 
contains no large logarithm. Renormalization group can be used to obtain C m 
at g <C m: 


CM = exp 


7 m(a) 

a s (m ) ^ y 


( 4 ) 


Two-loop anomalous dimension y m = = 7i + 72 (fj ) 2 H-and (3- 

function [3 = — \ = (3i^ + P 2 {^) 2 + • ■ ■ should be used together with 

one-loop terms in C rn (to) = 1 + C\ as ^ m ^ + • • ■ ||. Chromomagnetic interaction 
is the only term violating the heavy-quark spin symmetry |tJ] at the 1 /to level. 
Numerous applications of the lagrangian ([j] |j) are reviewed in (Q. 

In this paper, we obtain C m (g) at two loops from QCD/HQET matching. 
We consider scattering amplitude of an on-shell heavy quark with the initial 
momentum p\ = mv and a final momentum P 2 in a weak external field to the 
linear order in q = P 2 — pi- Similarly to Q, we use dimensional regularization 
(in d = 4 — 2e dimensions) and the background field formalism in which 
the combination gA l f \ is not renormalized. We consider QCD with rg massless 
flavours and a single heavy flavour Q ; the effect of a massive flavour with a differ¬ 
ent mass will be considered elsewhere. There are two possible effective theories: 
with loops of the heavy flavour Q and without such loops. Matching on-shell 
matrix elements produces finite results in both theories, in contrast to the case 
of a heavy-light current The HQET diagrams contain no scale and hence 
vanish in dimensional regularization, except for those with a massive quark 
loop (Section |2|). QCD on-shell diagrams were independently calculated with 
the REDUCE [0. [jj package RECURSOR (uj and a FORM (H) package @ 
(Section ||). Comparing the HQET and QCD matrix elements, we obtain the 
matching coefficient C m (m) and its anomalous dimension 7 (Section |]). 


2 HQET calculation 

The scattering matrix element of an on-shell quark in HQET up to 1/m level 
has the structure 

M = u v (q)gA“t a (V + C k Z ~+ C m Z^g [ ^>j u v { 0), (5) 


where qv = 0; e = Zq£ 0 , g = Zq/lq; £q and go are the bare proper vertex func¬ 
tions of the unrenormalized operators Q v ((vD) 2 — D 2 ) Q v and QyG^a^Qv. 
In HQET with m massless flavours in loops, all loop corrections vanish because 
they contain no scale: Zq = 1, £0 = 1, go = 1- 

In HQET with the heavy flavour Q in loops, the on-shell wave-function 
renormalization constant was found in Q : 


v ^ rj, 9o m ^(d- 1 )^- 2 )^- 6 ) 

1 — ^F-LF—r. —— 4n 


(47r) d 0 2(d — 5)(d — 7) 


( 6 ) 
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Figure 1: Diagrams for eo; + s means adding the mirror-symmetric diagram. 

Diagrams for £q are shown in Fig. [I]. We consider = 0, and therefore 
diagrams in which the background field is attached to the leading-order HQET 
vertex ig^v 11 do not contribute. Diagrams in which the background field is 
attached to a four-leg 1/m vertex do not depend on q and do not contribute to 
the structure linear in q. The diagrams Fig. jlje, h have zero colour factor. A 
method of calculation of scalar integrals in such diagrams was proposed in 
we average over directions of v in scalar integrals and obtain 


( vk) n 


r(j) 

r(i) r(^) 


( fc 2)n/2 


( 8 ) 


for even n (positive or negative) and 0 for odd n. After that, two-loop massive 
bubble integrals remain. An explicit formula for them can be found in (|l6t . We 
perform calculations in an arbitrary covariant gauge, and check that the sum 
of diagrams is gauge-invariant. Only the diagram Fig. |lja contains the colour 
structure CfTf ; it is easy to see that this contribution is exactly compensated 
by Zq (|j). In order to calculate epj we can contract the diagrams with q in 
the polarization index of the external gluon, and extract the q 2 part. In the 
backgroung field formalism, such external-gluon insertions produce the differ¬ 
ence of the propagators with the original momentum and the momentum shifted 
by q (after separating colour factors). Let’s consider the diagrams Fig. jljc, d, 
g first. Most of the terms in the sum cancel each other, leaving the difference 
of two terms: the diagram without the external gluon, and the same diagram 
with the heavy-quark momentum shifted by q. This shift does not influence the 
propagator; the kinetic energy vertex contains a term linear in q , which vanishes 
after the loop integrations. The sum of the diagrams Fig. 0b. f gives a similar 
difference of two diagrams without the external gluon. The difference of the 
two-leg kinetic energy vertices now contains a q 2 term, which is multiplied by 
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the same integral as in (^). Finally, we obtain 

£ = 1. (9) 

This fact is crucial for the proof of ©>■ Note that the “all-order” proof of 
the reparametrization invariance Q ignores all massive loops in HQET (of the 
external heavy flavour or any other massive quark) and therefore is valid up to 
the one-loop level only. 

In the case of Aco the diagrams Fig. [jjj, f do not exist. The colour factors 
of Fig. [lje, h no longer vanish. Again, the only CfTf contribution of Fig. jd|a is 
compensated by Zq and we obtain 


A = 1 + CaTf 


—4 £ 

do™ e j2 
(4 ir) d 0 


(d — 2)(d 2 — 9d + 16) 
4(d- 5)(d- 7) 


( 10 ) 


3 QCD calculation 


The scattering matrix element of an on-shell quark in a weak external field has 
the structure 


M = u(p 2 )gA“t a 



(pi +P 2 Y 
2 m 



u(pi) ■ 


( 11 ) 


Up to the linear terms in q it is determined by the total quark colour charge 
e = Zq£ o( 0) = 1 and chromomagnetic moment p = Zq/j, o( 0). The on-shell 
wave-function renormalization Zq was obtained in 0 - Let us see why e = 1. 
Diagrams for the bare proper vertex Aq in the constant background field A M can 
be obtained from the diagrams for the bare mass operator So (p) by shifting p: 
So (jp + A) = So(p) — Aq Afj,. The term £0 is separated by the 7 -matrix projector 
£0 = 1+\Tr Agt^^+l). On the mass shell, 1 Tr(S 0 (mu+A) — So(mt)))(^+l) = 
(1 ~Zq )vA. This gives £0 = Zq 1 . This argument is valid for both an abelian 
and a nonabelian background field A. 

We calculate bare proper vertices £ 0 , po on the renormalized mass shell. 
Therefore, it is convenient to use the pole mass m in the lagrangian, and to in¬ 
corporate the vertex produced by the mass counterterm Am fl(ij| . Diagrams for 
the proper vertex can be obtained from those for the mass operator by inserting 
the background field vertex in all possible ways (Fig. ||). Using integration by 
parts recurrence relations @• 0 . all two-loop on-shell integrals can be reduced 
to Iq and 


1 f d d kd d l 

n d J (-k 2 )(-(l - fc) 2 )(l - (v + l) 2 ) 

4(2d - 7) r 2 (i - £)r(i + 2£)r(i - 4 £ ) r2 
(d - 3)(3d - 8)(3d - 10) T (1 + £)r(l - 2 e )r(l - 3e) 1 ’ ’ 

1 f d d kd d l 

n d J (1 - (v + k) 2 )(l - (v + l) 2 )(l — (v + k + l) 2 ) 

3(d — 2) 2 (5d — 18) 2 n d — 4 16(d-4) 2 

2(d-3)(3d-8)(3d-10) 0 2d - 7 1 (3d - 8)(3d - 10) 1 ’ ’ 
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Figure 2: Diagrams for the QCD proper vertex. 
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where (T|, |l3 


/(e) = J +0(e), / = 7 r 2 log 2 — ^C(3)- 


(13) 


The diagrams Fig. |c with massless quarks (colour structures CfTftii and 
CaTftii) and Fig. Ed, e contain only I\. The integral I 2 is contained only in 
Fig. [Ic with the heavy quark (colour structures CfTf and CaTf ) and Fig. [|g. 
The diagram Fig. @h3 has zero colour factor. 

We perform all calculations in an arbitrary covariant gauge, and check that 
the d-dimensional results for e 0 and po are gauge-invariant. We check that £0 = 
Zq ; the same equality holds if we use colour factors for an abelian background 
fields in £ 0 . Moreover, this is true for each group of diagrams obtained from a 
single diagram for the mass operator: the sum of these diagrams, with each set 
of colour factors, gives the contribution of the original diagram to Zq 1 . This 
provides a strong check of our procedures. The programs for calculating /iq are 
obtained by replacing only the 7 -matrix projector and hence are reliable. If we 
use colour factors for an abelian background field, we reproduce the heavy-quark 
magnetic moment (llj; if the dynamical gauge field is also abelian, the classic 
result for the QED electron magnetic moment j^oj is reproduced (in dimensional 
regularization, it was discussed in |n]]). In view of all these checks, we are 
confident in our final result for the heavy-quark chromomagnetic moment. 

The full d-dimensional result is presented in the Appendix. Expanding it in 
£ and re-expressing it via a s (/z), we obtain 


A* = 


1 + T^('-5A) 


2Cf(1 + 4e) + Ca ( — + 2 + — —s + 2e 


-\Le 


C'f ( - 8 / + y 7T 2 - 31 ) + C F C A 


281 
3 £ 


■ z . 7 1 1011 4 t 3 

+ ci -^ + —- + - 1 -- 7T 


3 £ 


9 £ 

1 196 


1057 

~tT 


CfTfui — --— + CaTfUi — --s —-- -7T —— 


3 £ 


.CrTrl-* 1 —™*'*™ 
1 3e 3 9 


2 1 371 


371 


3 £ 2 9 £ 

4 1 


27 


605 

1 T 


(14) 


_ m , - * 8 1 8 2 370 

+ A F * 3£ 2 3£ + 9^ 27 


where L = log m//z, (3 1 = ^ Ca — nj = ni + 1. 


4 Results 

The coefficients Ck, C m in the HQET lagrangian ([|) are tuned in such a way 
that the full QCD matrix element (|ll]), expanded to linear terms in q , is equal to 
the HQET matrix element (||) . Infrared (or on-shell) singularities are the same 
in both theories. QCD and HQET spinors are related by Q u(pi) = it„(0), 
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u(p 2 ) = ( 1 + J- ) u v (q); the extra term with tfa should only be taken into 


account in zeroth order terms, but then it vanishes. 

We consider HQET with ni light flavours in loops first. Comparing the 
structure q we obtain Zjy 1 Ck = 1, and hence Ck(p) = 1 (§)• Comparing the 
structure [ 7 ^,$], we obtain Z^C m = p/p. We find Z m from the requirement 
that C m is finite. Terms 1/e 2 in it satisfy the consistency condition which is 
necessary for the anomalous dimension to be finite at e —■> 0. Terms 1/e give 
the anomalous dimension 


7m =2 C A ^ + ±C A (17C A - 13!>„,) (^) 2 . 
The chromomagnetic interaction coefficient at /i = m is 


(15) 


C m [rri) — 1 + 2(C F + Ca) 


a s (m) 

4n 


C 2 F ( -81 - 


20 


7r 2 - 31 


_ „ ,'4_ 4 2 269 

+ C F C A [- I + - tt + — 


( 4 17 SOS 

' 4 1 3 9 


27 


C F T F m ( 


rp , 4 2 299 

-C A T Fni (--n - — 


(16) 


+ C F T F [-™ir 2 + ™ 


+ C a T f n 2 - 


298 

"27" 


(~y 

\4tt / 


With N c = 3 colours, its numerical value is 

16 a s (rn ) 


C m (m) = 1 + + ( 21.79 - 1.91n,) (—) . 

3 7T V 7T / 


(17) 


The two-loop correction is large. The exact two-loop coefficient at ni = 4 is 
40% less than the expectation based on the naive nonabelianization i. e. it is 
not particularly accurate, but predicts the correct sign and order of magnitude. 
The heavy-quark loop contributes merely —0.10 to the bracket in (p~7|). 

If we now include Q-loops in HQET, we still have Ck{p) = 1 #• The 
anomalous dimension ( |l5| ) now contains rif = ni + 1 instead of ni . The chro¬ 
momagnetic interaction coefficient & has the coefficient of C A T F equal to 
^7r 2 — ^27, leading to 22.14 in the bracket in (|l7|). 

Our main results are the anomalous dimension ([Is]) and the chromomagnetic 


interaction coefficient at p = m (H3). If L = log m/p is not very large, the best 
approximation to C m (p) is the exact two-loop matching formula 

C m (p) = l+(Ci - 71L) s ^_ ^ + [Ca - (C171 + 72) L + 71 (71 - Pi) L 2 ] , 

V) 

in which all terms {a s / F ) 2 L 21,0 are taken into account, but {a s / F ) 3 L 3 and 
other higher order terms are dropped. Otherwise, it is better to sum leading 
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and subleading logarithms using 


C m (n) = 


u s {n) 

a a (m) 


- 7 l/( 2 / 3 i) r 


1 + Ci 


o-s (m) /3i72 - (3 271 (m) - (m) 


47T 


2(31 


47T 


( 19 ) 

These results can be applied to all cases of the spin symmetry violation, such as 
D-D* and B-B* splittings, 1 /to corrections in B —> D and B —> D* semilep- 
tonic decays, etc. [|). 

In the course of this work, we were informed by M. Neubert about the on¬ 
going calculation of the two-loop anomalous dimension of the chromomagnetic 
interaction operator by a completely different method 22 . Our result ( |l5| ) 
agrees with Q . 
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A Heavy—quark chromomagnetic moment 


The d-dimensional chromomagnetic moment has the form 

g 2 0 m - ld-2 ^ _ 4 ^ d _ _ ^2 _ 8d+ 14 j CA j I(j 


H — 1 + 


(47r ) d / 2 4 d — 3 


4 —4c 

i 9o m y at. 


(4tt) 




where the colour structures are Ci = Cp , C 2 = CfCa, C 3 = C'^, C 4 = CfTfUi , 
C 5 = CaTfUi, Cq = CfTf , C 7 = CaTf and 


do = 


Ji — 


(d — 2 )/q 


16(d - l)(d - 3) 2 (d - 4) 2 (d - 5)(d - 6)(d - 7) 

h 


J 2 = 


16(d-l)(d-3)(d-4)(2d-7)’ 16(d - l)(d - 4) 2 (d - 6) 

are chosen in such a way that the coefficients fly are polynomial in d. 
source of the following equation, presenting all nonvanishing coefficients a^-, 
was generated by a REDUCE program using the package RLFI by R. Liska 
(see 111, [Tlj| ): 


an = 4(d - 5)(d - 6)(d - 7)(d 8 - 32d 7 + 460d 6 - 3828d 5 + 19940d 4 - 66032d 3 
+ 135065d 2 - 155640d + 77356) 

a 12 = -16(2 d - 7)(2d 6 - 53d 5 + 557d 4 - 3005d 3 + 8828d 2 - 13450d + 8336) 
a 13 = 8(d - 6) (2d 6 - 60d 5 + 682d 4 - 3871d 3 + 11723d 2 - 18066d + 11120) 
a 21 = —2(d - 6)(d - 7)(2d 9 - 72d 8 + 1175d 7 - 11329d 6 + 70628d 5 - 293546d 4 
+ 809949d 3 - 1426799d 2 + 1454100d - 652948) 
a 22 = 4(16d 7 - 474d 6 + 5777d 5 - 37896d 4 + 145361d 3 - 327378d 2 + 402044d 
- 208160) 

a 23 = —4(d - 6)(5d 6 - 146d 5 + 1625d 4 - 9065d 3 + 27085d 2 - 41398d + 25424) 
a 31 = (d - 6)(d - 7)(d 9 - 35d 8 + 557d 7 - 5259d 6 + 32250d 5 - 132396d 4 
+ 362076d 3 - 633794d 2 + 642768d - 287288) 
a 32 = -(3d - 8)(6d 6 - 159d 5 + 1672d 4 - 9015d 3 + 26460d 2 - 40276d+ 24928) 
a 33 = 2(d - 3)(d - 6)(3d - 8)(d 4 - 23d 3 + 176d 2 - 550d + 596) 
a 42 = —64(d - 2)(d - 3)(d - 4) 2 (d 2 - 7d + 11) 
a 52 = 4(d - 3)(d - 4)(3d - 8)(2d 3 - 19d 2 + 57d - 56) 

061 = 4(d - 2)(d - 3)(d - 4)(d - 7)(3d 5 - 85d 4 + 875d 3 - 4123d 2 + 8898d 

- 7008) 

a 63 = —8(d - 4)(d 2 - 9d + 16)(d 3 - 13d 2 + 26d + 16) 

a n = -(d - 3)(d - 4)(9d 7 - 290d 6 + 3822d 5 - 26680d 4 + 106477d 3 

- 242974d 2 + 294012d - 145896) 

a 73 = 2(d - 4)(3d - 8)(d 4 - 16d 3 + 75d 2 - 96d - 28) 
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